function prepare_graph_combined(p_interest,ghat_MLE_tau25_unconstr,ghat_MLE_tau25_unconstr_noBerk,ghat_MLE_tau25_constr,ghat_MLE_tau25_constr_noBerk,ghat_MLE_tau50_unconstr,ghat_MLE_tau50_unconstr_noBerk,ghat_MLE_tau50_constr,ghat_MLE_tau50_constr_noBerk,ghat_MLE_tau75_unconstr,ghat_MLE_tau75_unconstr_noBerk,ghat_MLE_tau75_constr,ghat_MLE_tau75_constr_noBerk,data_directory, ref,pmin,pmax,y_interest,y_interest_univ)

	%==================================================================================================  
	fig_legend=figure;

		tau_text='50';
		eval(['plot(p_interest(1), ghat_MLE_tau' tau_text '_unconstr(1),  ''rv-'',''MarkerSize'',4,''MarkerFaceColor'',''r'',''MarkerEdgeColor'',''r'');']);
		hold on;
		eval(['plot(p_interest(1), ghat_MLE_tau' tau_text '_unconstr_noBerk(1), ''bo-'',''MarkerSize'',4,''MarkerFaceColor'',''b'',''MarkerEdgeColor'',''b'');']);
		eval(['plot(p_interest(1), ghat_MLE_tau' tau_text '_constr(1),  ''rs-'',''MarkerSize'',4,''MarkerFaceColor'',''r'',''MarkerEdgeColor'',''r'');']);
		eval(['plot(p_interest(1), ghat_MLE_tau' tau_text '_constr_noBerk(1), ''bd-'',''MarkerSize'',4,''MarkerFaceColor'',''b'',''MarkerEdgeColor'',''b'');']);
		leg_h = legend(strvcat('MLE, with Berkson', 'MLE, no Berkson','MLE (shape restricted), with Berkson', 'MLE (shape restricted), no Berkson'));
		legend_file(fig_legend, leg_h, [data_directory 'legend'], 'pdf');

	%==================================================================================================  
	fig_legend_bw=figure;

		tau_text='50';
		eval(['plot(p_interest(1), ghat_MLE_tau' tau_text '_unconstr(1),  ''kv-'',''MarkerSize'',4,''MarkerFaceColor'',''k'',''MarkerEdgeColor'',''k'',''LineWidth'',1.0);']);
		hold on;
		eval(['plot(p_interest(1), ghat_MLE_tau' tau_text '_unconstr_noBerk(1), ''ko-'',''MarkerSize'',4,''MarkerFaceColor'',''w'',''MarkerEdgeColor'',''k'',''LineWidth'',1.0);']);
		eval(['plot(p_interest(1), ghat_MLE_tau' tau_text '_constr(1),  ''ks-'',''MarkerSize'',4,''MarkerFaceColor'',''k'',''MarkerEdgeColor'',''k'',''LineWidth'',1.0);']);
		eval(['plot(p_interest(1), ghat_MLE_tau' tau_text '_constr_noBerk(1), ''kd-'',''MarkerSize'',4,''MarkerFaceColor'',''w'',''MarkerEdgeColor'',''k'',''LineWidth'',1.0);']);
		leg_h = legend(strvcat('MLE, with Berkson', 'MLE, no Berkson','MLE (shape restricted), with Berkson', 'MLE (shape restricted), no Berkson'));
		legend_file(fig_legend_bw, leg_h, [data_directory 'legend_BW'], 'pdf');

		%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
		% NOW COMBINE ACROSS THREE QUARTILES 
		%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
	 
	yy=2; % middle income group

	y_interest_sel=y_interest_univ(yy);
	sel=y_interest==y_interest_sel;

	tau_label_ydim=[6.95; 7.35; 7.70];

	fig_both=figure;
	hold on;

	tau_count=0;
	for tau=25:25:75

		tau_count=tau_count+1;
		tau_text=num2str(tau);
		yy_text=num2str(round(exp(y_interest_sel)));
		eval(['plot(p_interest(sel), ghat_MLE_tau' tau_text '_unconstr(sel),  ''ro-'',''MarkerSize'',4,''MarkerFaceColor'',''r'',''MarkerEdgeColor'',''r'');']);
		hold on;
		eval(['plot(p_interest(sel), ghat_MLE_tau' tau_text '_unconstr_noBerk(sel), ''bo-'',''MarkerSize'',4,''MarkerFaceColor'',''b'',''MarkerEdgeColor'',''b'');']);
		eval(['plot(p_interest(sel), ghat_MLE_tau' tau_text '_constr(sel),  ''rs-'',''MarkerSize'',4,''MarkerFaceColor'',''r'',''MarkerEdgeColor'',''r'');']);
		eval(['plot(p_interest(sel), ghat_MLE_tau' tau_text '_constr_noBerk(sel), ''bs-'',''MarkerSize'',4,''MarkerFaceColor'',''b'',''MarkerEdgeColor'',''b'');']);
		text(0.21,tau_label_ydim(tau_count),['(\tau = 0.' tau_text ')']);

	end % tau

	legend(strvcat('MLE, with Berkson', 'MLE, no Berkson','MLE (shape restricted), with Berkson', 'MLE (shape restricted), no Berkson'));
	grid on;
	xlabel('log price');
	ylabel('log demand');
	axis([pmin pmax 6.58 7.72]);
	name=[data_directory 'fig_combined_' ref '_yy' yy_text];
	savefig_fct_fullsheet(name,0);

	%-------------------------------------------------------------------------------------------------
	fig_both_BW=figure;
	hold on;

	tau_count=0;
	for tau=25:25:75

		tau_count=tau_count+1;
		tau_text=num2str(tau);
		yy_text=num2str(round(exp(y_interest_sel)));
		eval(['plot(p_interest(sel), ghat_MLE_tau' tau_text '_unconstr(sel),  ''ko-'',''MarkerSize'',4,''MarkerFaceColor'',''k'',''MarkerEdgeColor'',''k'',''LineWidth'',1.0);']);
		hold on;
		eval(['plot(p_interest(sel), ghat_MLE_tau' tau_text '_unconstr_noBerk(sel), ''ko-'',''MarkerSize'',4,''MarkerFaceColor'',''w'',''MarkerEdgeColor'',''k'',''LineWidth'',1.0);']);
		eval(['plot(p_interest(sel), ghat_MLE_tau' tau_text '_constr(sel),  ''ks-'',''MarkerSize'',4,''MarkerFaceColor'',''k'',''MarkerEdgeColor'',''k'',''LineWidth'',1.0);']);
		eval(['plot(p_interest(sel), ghat_MLE_tau' tau_text '_constr_noBerk(sel), ''ks-'',''MarkerSize'',4,''MarkerFaceColor'',''w'',''MarkerEdgeColor'',''k'',''LineWidth'',1.0);']);
		text(0.21,tau_label_ydim(tau_count),['(\tau = 0.' tau_text ')']);
	
	end % tau

	legend(strvcat('MLE, with Berkson', 'MLE, no Berkson','MLE (shape restricted), with Berkson', 'MLE (shape restricted), no Berkson'));
	grid on;
	xlabel('log price');
	ylabel('log demand');
	axis([pmin pmax 6.58 7.72]);
	name=[data_directory 'fig_combined_' ref '_yy' yy_text '_BW'];
	savefig_fct_fullsheet(name,0);
	%-------------------------------------------------------------------------------------------------
